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DIFFRACTION EFFECTS IN CERENKOV RADIATION 



John R. Neighbours and Fred R. Buskirk 
Physics Department 

Naval Postgraduate School 
Monterey, CA. 93940 

INTRODUCTION 

In previous work-'- Cerenkov radiation was considered for 
periodic bunches of electrons such as would be emitted by a 
typical travelling wave electron accelerator (Linac). Even in 
air, electrons from a 100 MeV Linac exceed the velocity of light 
in the medium, and Cerenkov radiation should be emitted at a small 
angle of less that 2 degrees. The radiation intensity was 
calculated in detail for microwave frequencies and some of the 
results were as follows: 

a. For bunches periodic in time with frequency v Q , 
radiation is emitted at v Q and harmonics therof, in contrast to 
the continuous frequency distribution observed for a single 
charge . 

b. If each bunch has a spatial distribution described by a 
charge density p ^ ( r ) , the radiated intensity is modified by 

the fourier transform of this charge distribution. 

c. At low frequencies such that the wave length of the 
emitted radiation is of the order of the bunch size, the electrons 
in the bunch radiate coherently. This leads to large enhancement 
factors for typical linac bunches consisting of 10® electrons, 

and in fact allows the radiation to be significant at microwave 
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frequencies. Destructive interference, described by the fourier 
transform of the charge density, decreases intensities with 
increasing frequency, until incoherent radiation takes over when 
the wave length of the radiation is much less than the electron 
spacing. 

d. If the emission region has finite length, which may be 
realized by passing the electron beam through a gas cell, the 
radiation propagation direction is not confined to a sharp 
Cerenkov angle q c , where cos e c is defined to be c/v, but is 
spread over a range of emission angles. This spreading effect 
depends only on the length of gas cell, and does not depend on the 
electron bunch structure or periodicity of the bunches. In this 
paper as well as 1, c is the speed of light in the medium. 

The spreading of the emission angle noted in (d) above, as 
calculated in the earlier paper, contained approximations which 
are invalid in many common situations, so that the situation is 
considered in more detail here. 



2 



GENERAL AND LIMITING RESULTS 



The power radiated from periodic bunches of electrons 
traveling at a velocity v greater that the speed of light c in a 
medium has been given in Equations 25, 26 and 27 of reference 1. 

For convenience they are written below using in slightly different 
terminology. The power radiated per unit solid angle dp/dft is 



T 00 




o 



where W(v,n), the power per unit solid angle radiated at the 
frequency v, was shown to be 



W ( v, n) = L 2 ^ 2 v q 2 sin29 I Pq(£) | 2 I 2 (u) 



The parameters describing the radiation are 
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where n x , n y , n z are components of the unit vector n in the 
emission direction, v is the frequency of the emitted radiation, 
L(=2Z') is the length of the gas cell, Pq(?) is the charge 
density distribution for one bunch, and the usual Cerenkov angle 
9 C is given by cos9 c = c /v. The bunch frequency V Q is 
equal to the electron velocity divided by Z, the electron bunch 
spacing ( v Q = v/z) and v the frequency of the emitted radiation 
will be a harmonic of v Q . 

The results above are general ones for the emission of 
radiation from periodic bunches of electrons. In (2) the factor 
L^i(u)^ is identical to the result obtained for the 
calculation of elementary Fraunhofer diffraction from a slit 
across which the phase varies linearly such as plane waves 
striking the slit at an angle to the normal. However in (2) the 
diffraction angle is not the usual one measured from the normal to 
the radiating line source. Here it is more convenient to measure 
the diffraction angle 9 from the electron beam line so that it is 
the complement of the usual one. 

The sin^9 factor arises from cross products used in 
calculating the Poynting vector. It is just the angular factor in 
the usual expression for the power radiated from a dipole oriented 
along the electron beam. Thus the expression (2) for W can be 
interpreted as the interference of dipole radiators whose phase 
varies linearly with position along the beam line. 

In the expression (2) for the radiated power per unit solid 
angle, the factor I(u)^ is strongly peaked at u = o. (When 
u = o, the radiation angle is equal to the Cerenkov angle 9 C .) 
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In an earlier paper W was integrated over solid angle to obtain 
the total radiated power, assuming that I was so strongly peaked 
that the other functions in the integrand, namely sin^B and 
p' 0 (k), could be evaluated at 9 = 9 C and taken outside the 
integral. Fig. 1 of the previous paper shows the behavior of I 2 
as a function of 9, and there is defined a width T = AG = ir/2Z'. 
If, in fact, the sin9 and Pg(k) factors are slowly changing 
through the peak in I, the radiation is concentrated near the 
Cerenkov angle 9 C and the total radiated power is given by Eq. 

(32) of the earlier paper, which is 

?, = ojo) v sin^9 I p ' (k) I (4) 

0) I 4 IT o c 1 ~ o ‘ 

where N is the number of pulses in the interacting region. 

When I (u) is not so strongly peaked, then the variation of 
the other terms sin^B and Pq( 1<) in the expression for W must 
be considered. In general both factors should effect the position 
and height of the maximum in W with the sin^B term being 
predominant at smaller angles. 

It is convenient to obtain a more specialized expression for 
W, and to that end, to introduce the harmonic number j = oo/a) 0 
and to assume particular a charge density distribution } 

, 2 . 2 ~ 2 , 2 _ 

i „ / b k 2. a k . 2 n , /c . 

p^(k) = q Exp ( — j cos 9 - — ^ sm 9) (5) 
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which corresponds to a spatial charge distribution defined by a 
gaussian 



x 



p o (r) = C 1 Ex P ( — 2 ■ “2 

a a 



(6) 



where is a normalization constant related to the total charge 
per bunch and a, b are size parameters for the radial and 
longitudinal dimensions of the bunch. 

Then (2) becomes 

W(v,0) = — — q 2 L 2 

2c < 7 ) 



- , - b 2 k 2 2 a 2 k 2 . 2 fl . 

x < j 2 I 2 (u) sin 2 0 Exp ( cos 0 _ sin 0), 



For a narrow electron beam, a, the parameter describing the radial 
dimension of the bunch can be neglected in (7), so that W can be 
written 



W ( v , 0 ) = Wj (0) = Q Dj (0 ) 



( 8 ) 



where the radiation function is 



Dj (0) 
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the constant Q is 



Q = 



4 2 

2yv o q 2 
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and the diffraction variable is 



u = — (cos0 c - COS0) (10) 

and k z has been written for k cos0. 

It is worthwhile noting that the frequency of radiation 
enters into W not only as the harmonic of the fundamental 

frequency but also through the wave vector k (in the medium) 
in the expression (10) for the diffraction variable. 
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MICROWAVE CERENKOV RADIATION IN A GAS 



Equation (5) is generally applicable for gaussian electron 
bunches moving in a dielectric at velocities greater than the 
speed of light and (8) is the special case for a fine beam. The 
central result of this work is a detailed discussion of the 
results predicted by (8) for the radiated microwave power in a 
gaseous medium. 

The Cerenkov angle defined by cos 9 C = c/v = 1/nS is small 
for most gases because n, the index of refraction of the gas is 
very close to one. Consequently the dependance of 9 C on the 
electron velocity is very slight since 6 must also be close to one 
in order to attain the Cerenkov condition. Thus even for 
extremely relativistic electrons traveling in a typical gas, 9 C 
is always small (less than 2 degrees). 

In many situations, the factor (sin u)2/ u 2 appearing in 
(8) and (9) will be essentially a delta function of u, and the 
radiation will appear at or very near the Cerenkov angle 9 C . 

This limiting behavior occurs as kL ^ “ , with a result given by 
(4) . 

Inspection of (10) shows that u = o when 0 = 9 C ; a 
condition which leads to a maximum in the diffraction function 
I^(u) = (sin u/u)2. Further inspection shows that if kL is 
large, then u will be large at angles significantly different from 
9 C , leading to I^(u) varying rapidly as a function of 9 . i n 
the limit as kL + “ the diffraction function behaves like a delta 
function, having only the central maximum at u = o. This limiting 
behaviour occurs regardless of the particular value of Q c . 
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For finite kL, the radiation is emitted with a finite range 
of angles centered about an angle other than 0 C . Although the 
proper variable for discussing the radiation pattern is the 
product kL/ it more convenient to discuss predicted results in 
terms of the variation with either k alone or only L alone. 

In the range of 0 in the vicinity of 0 C , the sin^0 factor 
in (8) changes very rapidly although p^Ck) changes very 
little. Consequently the radiation then becomes not only smeared 
about the Cerenkov angle, but the angles greater than 0 C account 
for more radiated power than the angles between 0 and 0 C . The 
result is that not only is the radiation distributed in angle 
about 0 C , but the distribution is distorted so the the peak of 
intensity may occur at several times 0 C - 

To calculate and plot expected results a set of parameters 
were chosen close to those available in an experiment. A 
dielectric constant of e = 1.000536 was assumed for air giving an 
index of refraction of n = 1.000268. The electron parameters were 
assumed to be those for a 100 MeV linac: fundamental frequency 

V Q = 2.85 GHz and a gaussian bunch parameter of b = 0.24 cm. 

The gas cell length was assumed to be 90 cm - also experimentally 
attainable . 

Under these conditions, the calculated behavior of D j (0 ) is 
shown in Fig. 1 for the two harmonics j = 3 and j = 5. The 
expected Cerenkov angle is calculated to be 0 c = 1.29° and the 
striking result is how far the maxima in the high intensity 
(first) lobes are displaced from 0 C . Each harmonic is displaced 
differently with the lower one displaced further from 0 C . In an 



9 



earlier paper the radiated power was evaluated assuming the 
sin^e factor to be constant; but even in that approximation it 
was noted that higher harmonics would be spread less from 9 C . 

Further consideration of Fig. 1 also suggests the large 
amount of total power radiated into larger angles, this latter 
effect being enhanced by an additional sin8 factor (from the solid 
angle) which multiplies (8) in the expression for total power. 

To illustrate, also, that the spreading of intensity about the 
Cerenkov angle is an interference effect associated with the 
finite gas cell length, Dj(0) is calculated for several cell 
lengths, 70, 90 and 150 cm (usually 90 cm is assumed elsewhere in 
this paper) for the harmonic j = 3 in Fig. 2. For longer cells, 
the peak in the main lobe moves toward 9 C but the approach is 
very slow. However, it is understood rather easily: the first 

null occurs when the Huygens waves emitted from the front and rear 
of the cell (of length L) differ in phase by 2 tt, which gives 
simply the condition (see appendix B) 

— « (cos9 - cos6) = 2 tt ^ 1;L ^ 

c c 



Then the reason for the large breadth of the lobe is that for the 
assumed condition, cos9 varies slowly so that a large change in 9 
is required for the 2 tt phase shift required by (11) for a null. 

To show what might be observed in an experiment. Fig 3 
displays the sum of Dj(0) for harmonics j = 3, 4, 5, and shows 

the expected washing out of diffraction zeros which occur at 



10 



different emission angles for different harmonics. These 
harmonics would be seen by an X-band detector, with sensitivity 
from 8 to 12 Ghz, assuming V Q = 2.85 Ghz for an S-band Linac. 

The effect of adding the three radiation functions is to enhance 
and shift the first lobe, and to smear out the interference 
effects at higher angles since the individual nulls in the Dj 
occur at different angles. 

As the product kL increases the position of the maximum in 
the radiation pattern decreases smoothly to approach the Cerenkov 
angle. This behaviour is shown in Fig 4 where L is increased 
keeping the frequency (and thus k) fixed at the value for the 
third harmonic. The approach to 0 C is smooth but significant 
differences are to be expected even at very long cell lengths of 
100 m. 
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POWER RADIATED 



In the earlier sections it was shown that power is radiated 
in a range about the Cerenkov angle, as is shown in Fig. 1 and 
Fig. 2. It is interesting to compare the total power radiated 
with the predictions of the approximation used in the earlier 
paper. In Fig. 5, two measures of power radiated are plotted as 
functions of frequency: (A) the total power and (B) the power 

integrated over the main diffraction lobe. Also shown is the peak 
power per unit solid angle in the main lobe. These display a rise 
at low frequency characteristic of Cerenkov radiation followed by 
a fall off at high frequencies, associated with the form factor of 
the bunch. The total power has a shape similar to that merely 
sketched in Reference 1. The peaks in (A) and (B) occur at 
different frequencies and the approximate expression (4) (Eq. 32 
of ref. 1) would peak at still a different frequency. 

A most surprising result is obtained when the total power and 
the power in the main lobe are compared to the approximate value. 
These results are given in Table I ; additionally, the ratio of the 
total power to that predicted by the approximation is tabulated in 
the last column. This ratio has values of about 50 for the lower 
harmonics, and gradually falls off at the higher frequencies. 
Noting from Fig. 5 that the power spectrum peaks at about the 
third or fourth harmonic, the above ratio has values from 35 to 
43, meaning that the total Cerenkov radiation per unit length, 
under assumed conditions would be larger by a factor of about 30 
compared to what is expected from an infinite gas cell. 

The physical reason for this large increase in the power 
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radiated from a finite cell or path length may be illuminated by 
examining (2). The factor Pg(k) = q F(k) describes the 
distribution of charge in the bunch and merely allows coherence at 
low frequencies, but partial cancellation as the wave length of 
emitted radiation approaches the bunch length, and does not 
contribute to the power increase. The latter effect is governed 
by the two factors sin^G and I^(u). I^(u) through its 

argument u governs phase matching of the electron and the emitted 
wave from the beginning to the end of the gas cell (of length L’). 
Perfect phase matching occurs for u = 0 or cos8 = cos9 c = c/v. 

But for a finite cell, perfect phase matching is not required; as 
9 changes from 9 C , the radiation from the various parts of the 
cell tends to cancel. In particular, a null should occur when the 
radiation from the front and rear elements are out of phase by 2t. 
This phase difference of 2^ is obtained at an angle 9 N for which 
kL = 2 tt, (see appendix 1) and I(u) = 0. The fact that u contains 
cos 9 and 9 is near zero for a gas results in a relatively broad 
peak. This breadth, combined with the sin^ 9 factor in the 
intensity results in the large intensity emitted at angles larger 
than 9 C , between 9 C and the angle for the first null, 9^. 

Also the sin^9 factor accomplishes the increase in power, even 
if 8 C is small enough so that the null for 9 less than 9 C 
(i.e. 2u = — 2 tt ) is out of the physical range cos 9 _< l. 
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CONCLUSIONS 



The object of this paper is to consider in more detail the 
phenomenon predicted earlier, that the Cerenkov radiation, in 
cases where the interaction region is finite, should be broadened 
about the Cerenkov cone. For the case of a gas, for which the 
Cerenkov angle is a few degrees, the broadening is asymmetric 
about the Cerenkov cone, and the radiated power displays various 
interference lobes. Considerable power (about half in the case of 
a gas) appears outside the main lobe, and the latter is peaked at 
angles much larger than the Cerenkov angle. 

An even more interesting result is that the Cerenkov power 
radiated per unit path length is increased by large factors 
approaching two orders of magnitude for a gas cell of finite 
length compared to one of infinite length. 

This increase in power is associated with the finite length 
of the radiating medium only. The effect should be approximately 
the same for either single or periodic electron bunches, and 
should occur whether the bunches are effectively point charges or 
have significant size. 

To recall results obtained earlier, periodic bunches produce 
Cerenkov radiation at harmonics of the bunch frequency, while a 
single bunch will radiate with a continuous frequency 
distribution. The effect of the bunch size is reflected in the 
factor of the fourier transform of the charge in the bunch, which 
causes the radiation to fall off at high frequencies while at low 
frequencies all charges in the bunch radiate in phase. 
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The results outlined previously might occur in other 
situations or applications. The distinct phenomena are: (A) 

Coherent radiation by an electron bunch for wavelengths larger 
than the bunch, (B) Radiation at harmonics of the bunch frequency, 
(C) Smearing of the Cerenkov angle and (D) Asymmetric smearing of 
the Cerenkov angle with an increase in the power radiated. C and 
D should occur for finite radiator length; D should occur only if 
0 C is very small. 

All of the above are easily realized in the microwave region, 
where appropriate bunching occurs for travelling wave 
accelerators. Possibly useful sub-millimeter wavelength radation 
could be realized. 

In the optical range, C and D could occur for radiating cells 
which could be as large as 100 y meter in length. A, B, C and D 
together could only be realized with small scale bunching that 
occurs in a free electron laser. 

For the X-ray region it is impossible to envision bunching 
fine enough to accomplish A and B. The dielectric constant is 
less than unity, so that Cerenkov radiation does not generally 
occur. It might be possible to accomplish C and D if small 
spectral regions could be found in a medium for which the Cerenkov 
condition is satisfied such as near an atomic resonance^. Then 
greatly enhanced radiation could be produced by using a thin 
layered structure for a medium. 
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Table I 



Total radiated power as determined by numerical integration and 
by the approximation of reference 1. Power in the first lobe of the 
radiation patttern is also included- Parameters are as described in 
the text. Units are arbitrary. 



ABC A/C 



Harmonic 


Total Power 

(Numerical 

Integration) 


Power in 1st 
Lobe (Numer- 
ical Inter- 
gration) 


Total Power 
(Approxima- 
tion) 




n = 1 


11.86 


6.65 


.1174 


101.02 


2 


13.53 


6. 59 


.228 


59.34 


3 


14.12 


6.34 


.324 


43.58 


4 


14.12 


5. 96 


.402 


35. 12 


5 


13.71 


5.49 


.459 


29.87 


6 


13.52 


5.01 


.491 


27.54 


7 


12.06 


4.37 


. 512 


23.55 


8 


10.99 


3. 78 


.492 


22.34 


9 


9.86 


3.20 


.464 


21.25 


10 


8. 72 


2.66 


.425 


20.52 



16 



REFERENCES 



1. "Cerenkov Radiation from Bunched Electron Beams", 

F. R. Buskirk and J. R. Neighbours, Naval Posgraduate School 
School Report Number NPS-61-83-003 (October 1982) 

2. "Cerenkov Radiation from Periodic Electron Bunches", 

F. R. Buskirk and J. R. Neighbours, to be published in 
Physical Review A . 

3. Bazzlev, V. A., V. I. Glebov, E. I. Denisov, N. K. Zhevago, 

M. A. Kumakhov, A. S. Khlebnikov and V. G. Tsinoev. Sov. Phys. 
JETP 54(5) , 884 (1981) . 



17 



ACKNOWLEDGEMENT 



This work was partially supported by the Naval Postgraduate School 
Foundation Research Program. 



18 



ft) 4 * 




ANGLE 

Figure 1 . Radiation function Dj for the harmonics j = 3 and 
j = 5 as a function of emission angle. Vertical scale is 
rbitrary. For j = 3 the first maximum occurs at 9.8° and the 
irst minimum occurs at 16.2°. For the 5th harmonic the values 
are 7.6° and 12.6° respectively. Fundamental frequency is 2.85 
GH Z . Electron beam parameters and index of refraction are as 
described in the text. The emitted power per unit solid angle is 
equal to a constant multiplying Dj as stated in (8). 
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0 10 20 
ANGLE 

Figure 2 . Third harmonic radiation function 03 ( 6 ) as a function 
of emission angle for gas (air) cell lengths of 70, 90 and 150 cm. 
Vertical scale is arbitrary, and the same for all three curves. 

The peaks occurs at 11.1 9.8 and 7.6 degrees respectively. 

Electron beam parameters as given in the text. 
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ANGLE 



Figure 3 . D s , the sum of Dj(©) for j = 3, 4, 5 as a function 

of©. Electron beam parameters and index of refraction are as 
described in the text. The emitted power per unit slid angle is 
eaual to a constant times D„ . The maximum in D_ occurs at 
©= 8.30°. 
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Fiaure 4. Dependance of the first maximum in D3(<9) as a 
function of gas (air) cell length. 



.• t 
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HARMONIC 

FIGURE 5 . Radiated power as a function of linac harmonic number. 
P is the total radiated power obtained by integrating Dj from 0 
to TF • Pi is the radiated power in the first lobe of the 
radiation pattern. D m is the maximum value of Dj,* i.e., the 
value of Dj at<9 m . All calculations are for L = 90 cm and 
other parameters as given in the text. The vertical scale is in 
cm^. Actual powers can be obtained by multiplying by 2 ttQ. 
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Appendix A 

Discussion of the Power Increase . 

An approximate expression for the power radiated is 
developed. It is the object of this section to explain the large 
increase in power found when the present results are compared to 
the results of the earlier paper. We start from the expression 
for the power radiated per unit solid angle, 

42.22 2 

2yv V rir 2 sin u 2 - 

W = 2 sin^9 — F^(k) (A1) 

c 4 u 



where 



u 



kL 

2 



(cos0 

c 



COS0 ) 



(A2) 



The total power radiated at a given frequency v, where 
v = jv Q , is 



P u = | / Wd « = 2u / Wd(cos0) 



(A3) 



- V 



where 



A 2 = 



. 4 2 _ 

4yv q 2 ^ r 
o 

c. kL 



(A4) 



and 



• 2 T 2 r . 2 

3 L sin u 2 ^ 

I = I Sin 0 — 2 — F (k) du 

u 



(A5) 
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Approximation 1. 

If we assume sin^ 0 is constant and equal to sin^ 9 C and 
assume F^(k) is constant, the integral of sin^u/u^ gives tt . 

Then the approximation to the integral, denoted by Ij_, is given 
by 



J 2 l 2 



I, = 



2 2 
F z (k) sin^9 



(A6) 



Approximation 2. 

A better approximation to A5 is obtained by still considering 
F^(k) constant, but considering sin^G to be a variable. Then 
we have for I 

. 2 t 2 r . 2 

D L - + j sin u 

I = F (k) sin 0 — = du (A7) 

4 J u 



where u is given by A2. 
Mow note that 



du 



kL 

2 



d (cos0 ) 



kL 

2 



sin6 



d0 



Thus 



. 2 t 2 _ T r .2 

D L 7 ^ kL I sin u 

I = F z (k) — / sin J 9 — = d0 

4 2 J u 



(A8) 



(A9) 
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When the integrand is plotted, it may be noted that the 
function to be integrated has a peak at Q m equal to about half 
the value of 0 n , where 0 n is the value of 0 at the first null 
beyond u = o. Furthermore sin^u/u^ i s found to be near unity 
at 0 m . 

We thus arrive at a slightly better approximation to I, 
denoted by l 2 > by letting sin 0 = 0 n /2, and the width of the 
peak be 0 n / 2 • 



kL Id 



. 2 2 

3 L 2 - 

I, = F^(k) 

2 4 2 \4 



& 



(A10) 



Then the ratio of ^2^1 




If we let 9 C = 1.29°, k correspond to the 3 r( ^ harmonic 

of 2.85 GH-, and L = 90 cm, we find 0 = 16.1°, and 

z n 



(All) 



= 19.7 



(A12) 



This value corresponds to the value 19.6 obtained by doing 
the integral. The significant result is that the large increase 
in power is associated with the broad diffraction lobe of this 
physical situation which couples with the sin^0 dependence of 
the Cerenkov power formula. Thus when the peak is shifted from 
1.29° expected for an infinite cell, large increases in power are 
expected. 
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Appendix B 

Angles for Minima in the Diffracted Cerenkov Radiation 
The radiation functon Dj(0) given by (9) and therefore the 
radiated power (8) go to zero when the diffraction variable u 
takes on successive values which are integer nultiples of ti . From 
(10) this condition is 



kL (cos0 c - cos0) = n 2tt 



(Bl) 



For a given frequency and cell length, Bl can be solved for the 
values of the angle which give a minimum in W. 

2 TT 

COS0 = cos0 c - n j^- ( B2 ) 



Or in terms of the fundamental frequency and the harmonic number 
of the Cerenkov radiation, the result is 

n c 

cos0 = cos© - — (B3) 

c j v L 
J o 

Figure B-l is a plot of the first four minimum angles for 
harmonics up to j = 20. Other parameters are as discussed 
previously. The figure shows the dimunition of the diffraction 
lobes as j increases (kL increases). 
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Figure 31 . Shift in angles at which the diffracted Cerenkov 
radiation is zero as a function of harmonic number of the 
radiation. Cell length is assumed to be 90 cm and otner 
parameters are as described in the text. 
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